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Abstract. We are conducting a systematic study of highly-ionised outflows in AGN 
using archival Suzaku data. To date we have analysed 59 observations of 45 AGN 
using a combined energy-intensity contour plot and Montecarlo method. We find that 
~ 36% (16/45) of sources analysed so far show largely unambigous signatures (i.e., 
Montecarlo proabilities of > 95%) of highly-ionised, high- velocity absorption troughs 
in their X-ray spectra. From XSTAR fitting we find that, overall, the properties o f the 
absorbers are very similar to those found recently by lTombesi et al.l (|2010ll201l|) with 
XMM-Newton for the same phenomenon. 



1. Introduction 

Recent systematic studies with XMM-Newton have found that as many as ~ 40% of 
local (i.e., z < 0.1) AGN show evi dence for highly-ioni sed, high- velocity material out- 



flowing from the central nucleus (|lbmbesi et alJ 12010). These outflows - which are 



often referred to as "disk-winds" owing to their origin on sub-parsec scales - are ob- 
servationally characterised by blueshifted Fe xxv He<2 and Fe xxvi Lya absorption lines 
and can have a kinetic energy flux which can exceed the 0.5% of Lhd threshold though t 
necessary to be significant in galaxy scale feedback scenarios dHopkins & ElvislEoiOh . 
Furthermore, it is often found that by integrating the kinetic energy flux over a realistic 
Super-Massive Black Hole (SMBH) e-folding time that the total energy output can be 
comparable to the binding energy of a typical 10 11 M^ galaxy buldge, which also sug 



gests that they may be energetically significant (e.g., iReeves etaDl2010l . iGofford etal 



l201lh . The study of these outflows is therefore key if we are to bridge the small- 
and large-scale processes occuring in galactic evolution and, ultimately, to help under- 
stand the apparent fundamental relationships o bserved between gala xy properties and 
the central nucleus (e.g., the Mbh ~ cr relation. iFerrarese et alJ [2006) 

To this end we are conducting a new systematic study of highly-ionised AGN out- 
flows using archival data for the Suzaku space observatory. The wide energy coverage 
(0.6 - 50.0 keV) offered by Suzaku, in addition to its high throughput and excellent 
energy resolution at FeK (~ 120 - 130 eV at 6keV), make it the ideal instrument 
to simultaneously disentangle the broad-band components of AGN spectra and hence 
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confidently search for any highly-ionised absorption lines. Here, we describe the pre- 
liminary results of this work. 

1.1. Sample selection 

The heterogenously selected Suzaku sample consists of all pointed AGN observations 
contained in the Data Archives and Transmission System (DARTS) hosted by the JAXA. 
Observations were further screened according to the following criteria: (1) Publically 
available as of September 2011; (2) a net exposure of > 50 ks to ensure sufficient pho- 
ton statistics in even intrinsically faint sources. Observations which were split into 
multiple sequences (e.g., variability monitoring campaigns) were included if the total 
net exposure was > 50 ks; (3) a minimum of 10 k counts in the 2 - lOkeV band to 
detect any highly-ionised absorption lines at a reasonable confidence level. 10k counts 
is roughly the minimum number required to detect absorption lines at Fe K Montecarlo 
significances of > 95%; (4) Compton-thin, i.e., observed Nu < cr? 1 ~ 1.5 x 10 24 cm" 2 , 
so that the intrinsic continuum emission, and any associated line-of-sight absorption, 
is directly observable. The resultant sample consists of sources which are predomi- 
nantly local but also includes PDS456 (z = 0.184) and gravitationally lensed quasar 
APM 08279+5255 (z = 3.9122). 



2. Analysis 

2.1. Fitting the continuum 

The X-ray spectrum of AGN typically consists of a primary power-law continuum, 
Compton reflection, and, in some cases, a soft excess. When absorbed by distant pho- 
toionised material these components can impart significant complexity into the Fe K 
band (e.g., the FeKa fluorescence line and associated edge). It is therefore impor- 
tant to first constrain the broadband continuum, and the effects of each of the spectral 
components to the Fe K band, before assessing for the presence of Fe xxv Hea and/or 
Fe xxvi Lya absorption. 

We utilise the full bandpass of the Suzaku XIS+PIN (i.e., 0.6 - 50keV) to con- 
strain the contribution of the various continuum co mponents, and, using table models 
generated using the XSTAR photoionisation code (Bautista &Kallman 2001), model 
the considerable effects that both fully- and partially-covering absorption can have on 
the broad-band spectral curvature and the resultant continuum parameters. In particu- 
lar, we give close attention to the FeK edge at ~ 7.1 keV and account for the effects 
that both Compton reflection and photionised absorption can have on the edge profile. 
The Fe Ka fluorescence line and associated reflection hump are characterised with RE- 
FLIONX whi ch calculates the r eflected emission from a face-on slab of Compton-thick 
material (iRoss & Fabianl 12005). In the soft X-ray band XSTAR absorption tables are 



added until an adequate fit to the spectrum is achieved. 
2.2. Searching for atomic features 

For a given best- fit model, i.e., accounting for all continuum components, soft X-ray 
emission lines and any warm absorption, we have searched for ionised atomic features 
in the Fe K band using unbiased and statistic ally driven energy-intensity contour plots 
similar to those used by (iTombesi et al.ll2010b . The contour plots are generated by step- 
ping a Gaussian in 25 eV intervals between 5.0 - 10.0 keV (rest-frame), with the line 
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normalisation allowed to adopt both positive and negative values. The resultant plots 
provide a powerful means to visually determine whether there are any possible atomic 
features present between 5 - lOkeV which warrent further invesitgation with Monte- 
carlo simulations (see section [231) . Energy-intensity contour plots were generated for 
each observation in the sample and all residuals with an F-test significance > 99% were 
then modelled with Gaussian profiles. If a spectrum appears to require a broadened 
Fe Ka component the residua were modelled with a broad Gaussian or a DISKLINE 
depending on which yielded the best phenomenological fit to the data. In these cases a 
second contour plot was generated - this time with the broad component included - to 
determine whether its presence had any effect on the absorption line detection. 

All residua which were found to have F-Test and Montecarlo significances of Pp > 
99% and Pmc > 95%, respectively, were self-consistently modelled using an XSTAR 
photoionisation table with a turbulent velocity, v tur b, roughly matching the full- width 
at half-maximum (FWHM) velocity width of the best-fit Gaussian line. The use of 
XSTAR enables us to measure the key parameters of the absorbing material, such as 
the ionisation parameter, column density, TVh, and outflow velocity, v ou t> relative to 
the observer. Importantly, this then allow us to indirectly probe how the likely energy 
budget of the outflow and assess whether it is likely to play a significant role in galaxy- 
wide feedback scenarios. 

In Figure Q] we show some examples of this fitting process when applied to two 
objects in the sample: Mrk766 and PDS456. The energy-intensity contours show 
that there is unambiguous evidence for blueshifted absorption at F-test confidences of 
> 99.99% in both sources. In Mrk766 (left panel) the two profiles are unresolved 
(p~ = 10 eV) and have an energy separation consistent with Hea and Lya. Both lines 
are fit with a single highly-ionised absorption zone described by Nu = (5. 6+ J 4) X 
10 22 cm" 2 and log£ = 3.67 ± 0.06, with a modest blueshifted velocity of v ou t = -(5.4 + 
0.9) x 10 3 km s" 1 relative to the systematic. Conversely, the absorption lines in PDS 456 
(right) are much more extreme. Our best-fit XSTAR model has two absorbers with 
tied N u = (1.50+g-^) x 10 23 cm" 2 and log£ = 4.26+°-} 3 , and outflow velocities of 
v ut,i = -(81.0 ± 2.1) x 10 3 km s" 1 and v out ,2 = -(90.0 ± 2.4) x 10 3 km s" 1 , respectively. 
Interestingly, this suggests that both lines are due to blends of Hea/Lya and, therefore, 
that the the disk-wind may be clumpy and have multiple components with different 
outflow velocities. PDS 456 is one of the most extreme objects in the sample. 

2.3. Montecarlo simulations 

Protassov et al.1 (l2002h showed that the F-test yields inaccurate statistics when dealing 
with complex spectral models containing narrow spectral lines. We, therefore, do not 
claim to robustly detect any lines based solely on the F-Test. Instead, we take the 
F-test probability, Pp» as a rough 'guide' when analysing the spectum and follow up 
all suspected absorption lines that have Pp > 99% with a Montecarlo (MC) analysis 
designed to determine their true veracity against random noise in the spectrum. 

The MC method is as follows: from the null model, i.e., no ionised lines present, 
we simulate T = 1000 fake XIS (FI) spectra using the FAKEIT command in XSPEC. 
In each of these simulated spectra a Gaussian profile with width equal to that of the 
suspected line is stepped every 25 eV between 5 - 9.5 keV (180 steps) with the centroid 
energy free to vary and the line normalisation allowed to adopt positive and negative 
values. A global A% 2 distribution was then generated by recording the maximum im- 
provement due to random noise, Ax* . , when stepping the line through each simulated 
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Figure 1. Data/model residuals and energy-intensity contour plots for Mrk 766 
(left) and PDS 456 (right). The data/model residuals are given to the best-fit con- 
tinuum model of each source with reflection modelled using the PEXRAV model to 
highlight the Fe Ka emission line. The continuous open contour corresponds to a 
Ax 2 = +0.5 and is intended to indicate the level of the continuum. From outer to in- 
ner, the closed contours correspond to F-test significances of 68%, 90%, 99%, 99.9% 
and 99.99%, respectively. The dashed vertical lines indicate the expected rest- frame 
energies of the Fe Ka, Fe xxv Hea and Fe xxvi Lya lines. 



pectrum. If a Gaussian profile in the real spectrum gives a A^ eal improvement, the MC 
probability of this line is then P M c - 1 _ (N/T), where N is the number of simulated 
spectra with Ay 2 . > Ay 2 , , and T is the total simulated spectra. 

r ^ noise ^real' r 



3. Results 

While work on the sample is ongoing we have fully analysed 45 AGN with this method. 
In these objects 16/45 (21/59 observations) show evidence for statistically significant 
(^mc > 95%) highly-ionised absorption lines. This corresponds to ~ 36% overall 
detect ion rate and is very similar to the global detection fraction found by lTombesi et al. 
(l2010h in local radio-quiet Seyfert galaxies. 

The observational signatures for the outflows falls into four categories: the most 
frequently observed absorbers are those consisting of either a single Fe xxvi Lya line 
or both a Fe xxv Hea/Fe xxvi Lya pair, which each accounting for 6/16 (~ 38%) of the 
detected outflows. A further 2/16 (~ 13%) are most likely due to Fexxv Hea blended 
with lower ionisation species of iron, and the final 2/16 have more than one absorption 
profile (i.e., PDS 456 in fig. Q]) which suggests the presence to a multi-velocity system. 
The current distributions of XSTAR absorber para meters are shown in fig. [2l O verall, 
the results are in good agreement with those of the lTombesi et al.1 (I2010l . l201lh XMM- 



Newton sample. The TVh values [panel (a)] in particular are very similar and, while 
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Figure 2. Histograms showing the distribution of absorber parameters: (a) 
column density in units of log(A^n/10 22 ); (b) ionisation parameter in units of 
log(£/ergcms -1 ); (c) outflow velocity in units of log(v ou t/c). In all panels the 
shaded /coloured area c orresp onds to this work while the dashed line shows the re- 
sults of Tomb esi et al.l {2010). 



the measured column densities cover a wide range (i.e., 10 21 < W H / cm" 2 < 10 24 ) 
there is a peak at (3 - 10) x 10 22 cm" 2 and a mean value of A^suzaku ~ 1 x 10 23 cm" 2 . 
This is consistent with the analagous values found with XMM-Newton. The distribution 
of ionisation parameter [panel (b)] covers the range range 3.1 < log^/ergcms" 1 ) < 
5.5 and, with a peak and mean value a t logQf /ergcms" 1 ) « 4.0, is again consistent 
with the results found by iTombesi et al.1 (l2010h . Despite this good agreement in terms 
of intrinsic absorber properties the current distrubtion of outflow velocities appears to 
differ somewhat, as is shown in panel (c). Whereas there is a sharp peak at ~ 0.1 c 
found with XMM-Newton, we find a smoother and more continuous range of velocities 
with Suzaku, ranging fro m 0.004 < v ou t/c < 0.5, with very few at the expected peak 
velocity of 0.1 c (e.g., see iKingll2010h . However, given that analysis of the sample is 
ongoing it is currently too soon to determine whether this is an intrinsic property of the 
absorbers or a result of, for example, low number statistics or intrinsic detector biases. 
The full Suzaku outflow sample and a thorough discussion regarding all of the work 
outlined here will be given in Gofford et al. (2012; in preparation). 
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